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O 1 I. INTRODUCTION 

CO 

Hadronic interactions are in principle explained by quantum chromodynamics (QCD). Such a first-principles de- 
i scription of the hadron-hadron interaction, however, is highly complicated, particularly at low energy, where QCD is a 
nonperturbative theory. In practice, therefore, effective hadronic Lagrangians are often used. The coupling constants 
at the hadronic vertices are then among the most fundamental quantities that should be computed from QCD. 
1 1 ' There is a long history of successful approaches of describing the two-baryon interaction using meson-exchange 
potentials. The values of the meson-baryon coupling constants have been empirically determined so as to reproduce 
the nucleon-nucleon (NN) p], Q , hyperon-nucleon (YN) [H and hyperon-hyperon (YY) interactions in terms of 
e.g. one-boson exchange (OBE) models. The scalar mesons play significant roles in such phenomenological potential 
models. In early OBE models for the NN interaction the exchange of an isoscalar-scalar "cr" meson with a mass 
V^j ' of about 500 MeV was needed to obtain enough medium-range attraction and a sufficiently strongspin-orbit force. 
It was only later understood that the exchange of a broad isoscalar-scalar meson, the e(760) QQ, simulates the 
exchange of such a low-mass "a" @|. 

Existence of the "cr" meson is expected also from chiral symmetry of QCD. The a meson appears as a chiral partner 
of the Nambu-Goldstone boson, w, and thus plays the role of the "Higgs boson" in chiral symmetry breaking. Recent 
analyses of ir — tt scattering have revealed a broad resonance at the mass around 600 MeV, called /o(600) by the 
Particle Data Group @, which is considered to play the role of e(760) in the boson exchange potential of the baryonic 
interactions. 

In terms of OBE models for two-baryon interactions, the dominant contribution to A A interaction comes from the 
scalar e(760) meson exchange ^3^3- The recent identification of \^He and the measurement of the AA pair suggest 
that the binding energy of AA (ABaa — TO MeV) is considerably smaller than the binding energy of NN 12] . This 
is in contrast to the outcome of the earlier measurement which is ABaa — 4.7 MeV [l^. This issue has also been 
examined within the framework of Nijmegen OBE potential model D (NHC-D) [l4|. In this model, the pseudoscalar 
octet 7T, 77, 7/ , K, the vector octet p, </>, u>, K* and the scalar singlet e are the exchanged mesons and the coupling 
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constants are fitted to data while other physical properties of the particles are taken from experiment. The estimated 
value of ASaa hi this model implies a rather strong attractive AA. However, in the Nijmegen soft-core (NSC) potential 
models 0,01, where there is a scalar nonet instead of a scalar singlet, we have much weaker attractive potentials than 
in the case of NHC-D in the AA systems. In this framework, since the scalar £(760) exchange plays the most crucial 
role in A A interactions, it is necessary to determine the AAe coupling constant in a model independent way in order 
to understand the role of e exchange in the strangeness S — —2 sector. 

The structure and even the status of the scalar mesons, however, have always been controversial [TEl IT(|. In the 
quark model, the simplest assumption for the structure of the scalar mesons is the 3 Po QQ states. In this case, the 
scalar mesons might form a complete nonet of dressed qq states, resulting from e.g. the coupling of the P-wave qq 
states to meson-meson channels (Tj 1 ]. Explicitly, the unitary singlet and octet states, denoted respectively by E\ and 
£g, read 

£i = {uu + dd + ss)/vo , 

£ 8 = (uu + dd- 2ss)/VQ. (1) 
The physical states are mixtures of the pure SU (3)-fl.avor states, and are written as 

£ = cos 9 S £i + sin 9 a e% , 

fo = - sin S £i + cos 9 S £ 8 . (2) 
For ideal mixing holds that tan# s = l/v2 or 6 S ~ 35.3°, and thus one would identify 

£(760) = {uu + dd)/V2, 

Jo (980) = -as. (3) 

The isotriplet member of the octet is aj'°(980), where 

a°(980) = (uu-dd)/V2. (4) 

An alternative and arguably more natural explanation for the masses and decay properties of the lightest scalar 
mesons is to regard these as cryptoexotic q 2 q 2 states In the MIT bag model, the scalar qq states are predicted 
around 1250 MeV, while the attractive color-magnetic force results in a low-lying nonet of scalar q 2 q 2 mesons 0,^3- 
This nonet contains a nearly degenerate set of I = and 1=1 states, which are identified as the / (980) and a* ,0 (980) 
at the KK threshold, where 

a°(980) = {sdsd - susu) / V2 , 

/o(980) = {sdsd + susu)/V2, (5) 

with the ideal-mixing angle tan# s = — V2 or 6 S ~ —54.8° in this case. The light isoscalar member of the nonet is 

£(760) = udud. (6) 

The nonet is completed by the strange member k(880), which like the £(760) is difficult to detect because it is hidden 
under the strong signal from the if* (892) [T1[T^|. We shall use in this paper the nomenclature (a '°, /o, c, n) for the 
scalar-meson nonet, where one should identify a = £(760). 

One way to make progress with the scalar mesons is to study their role in the various two-baryon reactions (NN, 
YN, YY). Our aim in this paper is to calculate the AAcr, SSa and EEcr coupling constants, using the QCD Sum Rules 
(QCDSR) method. The QCDSR method [2(j is a powerful tool to extract qualitative and quantitative information 
about hadron properties |2ll l22T |. In this framework, one starts with a correlation function, which is constructed 
in terms of hadron interpolating fields. On the theoretical side, the correlation function is calculated using the 
Operator Product Expansion (OPE) in the Euclidian region. This correlation function is matched with an Ansatz 
which is introduced from the hadronic degrees of freedom on the phenomenological side, and this matching provides a 
determination of the hadronic parameters like baryon masses, magnetic moments, coupling constants of hadrons and 
so on. 

There are different approaches in constructing the QCDSR (see e.g. Ref. for a review). One usually starts 
with the vacuum-to-vacuum matrix element of the correlation function that is constructed with the interpolating 
fields of two baryons and one meson. However, this three-point function method has as a major drawback that 
at low momentum transfer the OPE fails. Moreover, when the momentum of the meson is large, it is plagued by 
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problems with higher resonance contamination |23(. The other method that is free from the above problems is the 
external-field method. There are two formulations that can be used to construct the external-field sum rules: In the 
vacuum-to-meson method, one starts with a vacuum-to-meson transition matrix element of the baryon interpolating 
fields, where some other transition matrix elements should be evaluated [2]]]. (This is also the starting point of the 
light-cone QCDSR method.) In Ref. [24|, pion-nucleon coupling constant was calculated in the soft meson limit using 
this approach. Later it was pointed out that the sum rule for pion-nucleon coupling in the soft-meson limit can 
be reduced to the sum rule for the nucleon mass by a chiral rotation so the coupling was calculated again with a 
finite meson momentum |25| . These calculations were improved considering the coupling schemes at different Dirac 
structures and beyond the chiral limit contributions [26l [2?l l2p . In this paper, we calculate the baryon-sigma meson 
coupling constants, using the external field QCDSR method 29]. We evaluate the vacuum to vacuum transition matrix 
element of two baryon interpolating fields in an external sigma field and construct the sum rules. This method has 
been used to determine the magnetic moments of baryons [23. 130L I31L 13 ^ ■ the nucleon axial coupling constant [32LI33]] . 
the nucleon sigma term |34| . and baryon isospin mass splittings |35|. It has also been shown that at low momentum 
transfer, this method is very successful in evaluating the hadronic coupling constants. Recently, the NNa coupling 
constant, gNNa, was calculated using this method [3jj. It has also been applied, previously, to the calculations of the 
strong and weak parity violatin g pi on-nucleon coupling constants |37l I38L [39| and the coupling constants of the vector 
mesons p and ui to the nucleon [40|. 

In the SU(3) flavor symmetric one can classify the meson-baryon coupling constants in terms of two parameters, 
the NNa coupling constant, gMNa an d the Fj (F + D) ratio of the scalar octet, a s ,41] : 

1 2 

9NNa =g, g N Ne s = {±ot s - I) , g AAes =--j=g(l-a s ), (7) 

1 2 

gsB.es = --/=fl , ( 1 + 2a a ), sr S £ £8 = -j=g(l-a s ), g E Ba = g (2a s - 1) , 

2 1 

ffESao = 2 5 3AEa = # C 1 ~ a «) ' 9ANk = - g (1 + 2a s ) , 

guNK. = g (l - 2a s ) . 

Considering the mixing between the singlet and the octet members of the scalar nonet, one obtains for gBBo and 

9BBf , 

gBBa = cos 9 s gi + sin O s g B Be 8 > 

gBBfo = - sin 9 s gi + cos O s gBBs 8 , (8) 

where g\ = gBBa is the flavor singlet coupling, and 9 S is the scalar mixing angle. 

We shall first consider the sum rules in the SU (3) flavor symmetric limit to see if the predicted values for the meson 
baryon coupling constants from the sum rules are consistent with the SU(3) relations. We show that this is indeed 
the case which leads to a determination of the F/(F + D) ratio of the scalar octet. Furthermore, keeping track of 
these coupling constants with the SU{3) relations, we obtain the values of the other scalar meson-baryon coupling 
constants. For this purpose, we assume ideal mixing and make the analysis in both qq and q 2 q 2 pictures for the scalar 
mesons. As we move from the S — to the S = — 1 and S = —2 sectors, the flavor SU(3) breaking occurs as a result 
of the s-quark mass and the physical masses of the baryons and mesons. We also consider the SU (3) breaking effects 
for the sum rules to estimate the amount of breaking, individually for each coupling. 

We have organized our paper as follows: in Section^ we present the formulation of QCDSR with an external scalar 
field and construct the sum rules for the AAcr, SSer and EEcr coupling constants. We give the numerical analysis and 
discuss the results in Section ITTll Finally, we arrive at our conclusions in Section llVl 

II. BARYON SUM RULES IN AN EXTERNAL SIGMA FIELD 

A. Construction of the Sum Rules 

In the external-field QCDSR method one starts with the correlation function of the baryon interpolating fields in 
the presence of an external constant isoscalar-scalar field a, defined by the following: 

n B *(q) =i fd 4 x e l "- x (o\T[r, B (x)fj B (0)]\o^ , (9) 
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where tjb are the baryon interpolating fields which are chosen as |2l| 

Vs = ea6c[(sfC7 /i s fc )7 5 7 M u c ] , (10) 

r) A = (2/3) 1 ^h abc [(u^C^s b h 5 ^d c -(d r [C^ Sb h 5 ^u c } . 

for S, S and A, respectively. Here a,b,c denote the color indices, and T and C denote transposition and charge 
conjugation, respectively. For the interpolating field of each octet baryon, there are two independent local operators, 
but the ones in Eq. (|1C)|> are the optimum choices for the lowest-lying positive parity baryons (see e.g. Ref |42j for a 
discussion on negative-parity baryons in QCDSR). 

The external sigma field contributes the correlation function in Eq. @ in two ways: first, it directly couples the 
quark field in the baryon current and second, it modifies the condensates by polarizing the QCD vacuum. In the 
presence of the external scalar field there are no correlators that break the Lorentz invariance; however, the correlators 
already existing in the vacuum are modified by the external field: 

(qq)a = (qq) + g° q x<r{qq) , (u) 

(g c q(T ■ Gq) a = {g c qcr ■ Gq) + gqXG<^(g c qo- ■ Gq) , 

where only the responses linear in the external-field are taken into account. Here, g q is the quark-cr coupling constant 
and x and \G are the susceptibilities corresponding to quark and quark-gluon mixed condensates, respectively. In 
Eq. l(TT|l . (qq) represents either (uu) or (dd), as we have assumed that (uu) ~ (dd) and the responses of the up and 
the down quarks to the external isoscalar field are the same. Note that, here we assume ideal mixing in the scalar 
sector, that is, we take the sigma meson without a strange-quark content. Therefore, the sigma meson couples only 
to the u- or the d-quark in the baryon, where we take g^ = g% and g a s = 0. 

In the Euclidian region, the OPE of the product of two interpolating fields can be written as follows: 

n B!T {q)=Y d C°{q)O n , (12) 

n 

where C^{q) are the Wilson coefficients and O n are the local operators in terms of quarks and gluons. In order 
to calculate the Wilson coefficients, we need the quark propagator in the presence of the external sigma field. In 
coordinate space the full quark propagator takes the form: 

S q (x) = SW(x)+SW(x), (13) 

where, 



and 



i S^ a 



' - (0|T[ (Z a (x)g b (0)|0) o 

„• :ab ■ \n i sab 



S ab x 2 

192 



„q(T ■ Gq) 



(14) 



-^KbdcG^a^ ln(-x 2 ) - 29x3 J m q x 2 ln(-x 2 ) 



Att 2 x 2 32tt 2 



mq (qq)x + 1 ™ qS * (g c q<r ■ Gq)x + 0(a 2 , m 2 q ) 



48 



2 7 x 3 2 



i Si a)ab 



(0\T[q a (x)q b (0)\0) a 
S ab 



(15) 



9 a O 



ig ab S a b)( i S ab x 2 

\qq) x — 77r\qq) + wf — 7^{9cqa- ■ Gq)x 



48 

5 ab x 2 
192 



12 ' 2 7 x 3 2 
X G{g c qcr ■ Gq}] + Q(a 2 ) . 
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Here, G^ v is the gluon tensor and g 2 — Aira s is the quark-gluon coupling constant squared. Note that, in the quark 
propagator above, we have included the terms that are proportional to the quark masses, m qi since these terms give 
non-negligible contributions to the final result as far as the strange quark mass is considered. 

Using the quark propagator in Eq. one can compute the correlation function T\b<j(si)- The Lorentz covariance 
and parity implies the following form for Hbo-(<i)- 

n Ba (q) = (H^ + n« q) + (H^ + q)a + 0(a 2 ) , (16) 

where q = q^j^ is the four momentum of the baryon. Here Hg and Tl q BQ represent the invariant functions in the 
vicinity of the external field, which can be used to construct the mass sum rules for the relevant baryons, and Hg a 
and n?, denote the invariant functions in the presence of the external field. Using these invariant functions, one 
can derive the sum rules at the structures 1 and q. In Ref. [3(| it was found that the sum rule at the structure q 
for the NNa coupling constant is more stable than the other sum rule at the structure 1, with respect to variations 
in the Borel mass. Motivated with this, we here present only the sum rules at the structure q and use these for the 
determination of the coupling constants. 



B. A Sum Rules and AAa Coupling Constant 



We shall first present the sum rules calculations for AAc coupling constant in detail and in the next subsection, we 
shall give the sum rules for SScr and SScr couplings. At the quark level, we have for A: 



T[ VA (x)fj A (0)] 



2i 



^abc ^a' b' c' 



(Tr{Sf(xh u C[S b s b '(x)} T C 1 ,} l5 rS c /(xh^ 5 
+Tr{S c /(x) llJ C[Sf(x)} T C^h 5 rS a u a '(x) 1 ^ 5 

-757^r' (xh»C[Sf (x)] T C^Sf (xh"j 5 ) 



(17) 

Using the quark propagator in Eq. I|13[l . the invariant function at the structure q in the presence of the external 
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Act' 



field, n A is calculated as 



KM = 9° 



i 



(27T) 



4 4 777 - 

-a, (1 - /) ln(-q 2 ) - a*(l + 2/) - ( X + Xg)^ a 2 q (l + 2f) 
5 2 2 



12g 2 
1 

Qq 1 



(5m s - 3m g ) + 8(m s - m q )f 



(18) 



where we have defined / = fel — 1, a q = — (2n) 2 (qq) and {g c qc • Gq) = m 2 (qq). 

In order to construct the hadronic side, we saturate the correlator in Eq.@ with A states and write, 

IT (n\ <°I^I A ) /M„AS < A I^I°) 



(19) 



A 1 1V1 A 

where Ma is the mass of the A. The matrix element of the current ?7a between the vacuum and the A state is defined 
as, 



|?7a|A) = X A v , 



(20) 



where Aa is the overlap amplitude and v is the Dirac spinor for the A, which is normalized as vv — 2M A . Inserting 
Eq. i|20|) into Eq. (|19f) and making use of the isoscalar scalar meson-baryon interaction Lagrangian density 



we obtain the hadronic part as: 



£ = 9 A Act VV a , 



.. |2 q + M A q + M A 



q 2 -M 2 - 



q 2 -M 2 



(21) 



(22) 
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We have also contributions coming from the excitations to higher A states which are written as, 



AaAa« 



q + M A 
q 2 -M 2 K 



<?AA* 



g + M A * 
q 2 - M%. 



(23) 



and the ones coming from the intermediate states due to er-A scattering i.e. the continuum contributions. Note that 
the term that corresponds to the excitations to higher A states also has a pole at the A mass, but a single pole instead 
of a double one like in Eq. Ij22(l . This single pole term is not damped after the Borel transformation and should be 
included in the calculations. There is another contribution that comes from the response of the continuum to the 
external field, which is given as: 



-Ago b{s) 

2 — d(s - s )ds . 



(24) 



where sq is the continuum threshold, Asq is the response of the continuum threshold to the external field and b(s) is 
a function that is calculated from OPE. When Aso is large, this term should also be included in the hadronic part 

m 

Matching the OPE side with the hadronic side and applying the Borel transformation, the sum rule for AAtr 
coupling at the structure q is obtained as: 



A A M 2 

■| M 4 a q (1 - /) E£ + a 2 q (l + 2/) 



12 



M 2 mla q L- u ' 27 



(25) 



( X + XG)^a 2 (l + 2f)L- 2 / 2 ' + ±(2r 



18 



m , 



,)M 6 E^L- S / 9 



a q (2m s - m q ) M 4 £ A L^ 9 + ± [(5m. - 3m g ) + 8(r 



m q )f 



- 2 M A - AP 

"A A — — ffAAo- + OK 



A4 



9 a q 2<£ 



{s%y-Am s fa q 



M 2 L- 4 / 9 e ( M A- s o)/^ 2 



where we have defined A A = 327r 4 A A and M is the Borel mass. The continuum contributions are included by the 
factors 



E? 



e - s «/ M , 



(26) 



where s A is the continuum threshold. In the sum rule above, we have included the single pole contribution with 
the factors B\. The third term on the right hand side (RHS) of Eq. (|25|l denotes the contribution that is ex- 
plained in Eq. I|24|). Note that this term is suppressed by the factor e~( s »~ M A)/ M2 as compared to the single 
pole term. We have incorporated the effects of the anomalous dimensions of various operators through the factor 
L = \n(M 2 / Aq CD ) / \n(n 2 / Aq CD ) , where jU is the renormalization scale and Aqcd is the QCD scale parameter. 



and E Sum Rules and zJEo and EEcr Coupling Constants 



One can apply the method explained in the previous subsection for the S2<r and EEcr couplings and derive the 
corresponding sum rules. Using the interpolating fields in Eq. i|10fl . we obtain 



0) = 2ie abc e a ' b ' c ' 



Tr{ST {x) lv C[Sf (z)] T C 7m }757^ c Wis , 



(27) 



abc^a b c 



2ie aoc e 



Tr{Sf (xh„C[S b u b '( X )fC^} l5 rS? (x)Yj 5 . 



at the quark level for 5 and S, respectively. Using the quark propagator in Eq. I|13|) . the invariant functions at the 
structure q are calculated as: 



nL(g)=.< 



i 



^a q + a q ln{-q ) 



(28) 
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(2ny 



a i - ^9 + X)-^i a q - — a q - 2m q q ln(-g ) 
2m 

-2m«j x a 9 ln(-g 2 ) + a 2 



(29) 



The sum rules are obtained by matching the OPE side with the hadronic side and applying the Borel transformation. 
As a result of this operation, we obtain: 



^a, M 2 L- 14 / 27 - a q M 4 e M ^ M 
6 . 

= _Al ^ gu^ + — + ^ As= M 2 L -4/9 e (M=- s -)/M= 
" 9° q " .9? 2 5 g 



(30) 



and 



^ a 9 M 2 L- 44 ' 27 - (x G + X) ^ «^- 2/27 
-2 x m 9 a q M 4 E% L~ 4 / 9 + ~m q a 2 



2 M 2 L 4 / 9 +2m g M 6 £f£- 8 / 9 



3 Ml/M 2 



(31) 



-A 



B S ^ + |f[(^) 2 - 2m.(/ + l)a 9 ] M 2 i -4/9 e (Ml-^)/M 2 

5 g ^9q 



for SSer and ££cr couplings, respectively. 

We would like to note that the sum rule for SStr coupling constant at the structure q is independent of the 
susceptibilities x and Xg- Another feature of the sum rules above is that up to the dimension we consider, the terms 
involving the s-quark mass do not contribute to the OPE side. The contributions that come from the excited baryon 
states and the response of the continuum threshold are taken into account by the second and the third terms on the 
right-hand side (RHS) of the sum rules, respectively. 

For the sake of completeness, here we also give the sum rule for NNa coupling constant at the structure q [3(j: 



+2 m q M 6 E? L- 8 / 9 + 2 X m q a q M 4 E™ L" 4 / 9 + ~m q a 2 



G 

M 2 N /M 2 



2/27 



-M 4 a q E- + - 3 x M 2 a\ L 4 ' 9 - '-f M 2 a q L~ 14 / 27 ( X + Xg ) ^a 2 q L 

2 
3 

= -A?v ^£ 9nn<j + B N ^+ ( 4^ ^0 M 2 L -4/9 e (M^-^)/M= m 

9q 9q ^9q 

which follows from a choice of the interpolating field 0] 

Vn = e a b c [UaCj IJ ,Ub]j5l' J 'dc ■ (33) 

Note that, in Eq. I)32|) . the NNa sum rule in Ref. [3j| has been improved including the quark mass terms. 

Comparing the left-hand sides (LHS) of the sum rules in Eq. I|3U|I . Eq. (|31ll and Eq. (|32|l one can derive a basic 
relation between the SSer, SEtr and NNa coupling constants in the SU(3) limit, which is 



9nno 



(34) 



This relation is quite natural because for the SScr coupling only the it-quark propagator outside of the trace and for 
the EEcr coupling the it-quark propagators inside the trace involve the terms that are proportional to the external 
field. For the NNa coupling all the three quark propagators involve such terms and this implies that in the SU (3) 
and isospin symmetric limit, the relation in Eq. (|34|l holds. It is interesting to note that this relation can also be 
derived from Eq. Q and Eq. (jSJ) assuming the ideal mixing for the qq picture where the NNfo coupling vanishes. 



III. ANALYSIS OF THE SUM RULES AND DISCUSSION 



In this section we analyze the sum rules derived in the previous section in order to determine the values of the AAc, 
ES<7 and EEcr coupling constants. To proceed to the numerical analysis, we arrange the RHS of the sum rules in the 
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form 

/(M 2 ) = A B + B B M 2 + C B M 2 L- 4 /V M S- s *)/ M2 , (35) 
and fit the LHS to f(M ). Here we have defined 

A X2 M B 

A B = ~9BBcr ; 

9q 

Bb = ^, (36) 



.9, 



together with 
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C A = ^[(4) 2 -4m s / a ,], (37) 



C E = ^i[( So s ) 2 -4m s (/ + lK], 

for AAcr, SScr and EEcr sum rules, respectively. 

For the vacuum parameters, we adopt a q = 0.51 ±0.03 GeV 3 , and m§ = 0.8 GeV 2 [45|. We take the renormalization 
scale /i = 0.5 GeV and the QCD scale parameter Aqcd — 0.1 GeV. The value of the susceptibility \ nas been 
calculated in Ref. [2f| a s X = — 10 ± 1 GeV -1 . The value of the susceptibility xg is less certain. Therefore, we 
consider xg to change in a wider range. 



A. Scalar Meson-Baryon Coupling Constants in the 5(7(3) Symmetric Limit 

We shall first consider the sum rules in the SU(3) flavor symmetric limit, where we take m q = m s = and / = 0. In 
this limit we also set the physical parameters of all the baryons equal to the ones of the nucleon; Mb = Mjv = 0.94 GeV, 
As = \ N = 2.1 GeV 6 Hg, s$ = s$ . 

In Figs. 11131 we present the Borel mass dependence of the LHS and the RHS of the sum rules for AAcr, EEcr and 
SScr, respectively, for Sq = 2.3 and xg = X = ~ 10 GeV -1 . As stressed above, in the SU(3) limit we choose the 
Borel window 0.8 GeV 2 < M 2 < 1.4 GeV 2 which is commonly identified as the fiducial region for the nucleon mass 
sum rules 29] . It is seen from these figures that the LHS curves (solid) overlie the RHS curves (dashed) . In order to 
estimate the contributions that come from the excited baryon states and the responses of the continuum threshold, we 
plot each term on the RHS individually. We observe that the single-pole terms (dotted) give very small contributions 
to the sum rules except to the one for SScr coupling. The responses of the continuum thresholds (dot-dashed) for all 
the couplings are quite sizable. 

In order to see the sensitivity of the coupling constants on the continuum threshold and the susceptibility x, we plot 
in Figs. ^ and [S] the dependence of ghho/ 9q and gss^/gq on x for three different values of the continuum thresholds, 
Sq = 2.0, 2.3, and 2.5 GeV 2 , and taking x = Xg- One sees that these coupling constants change by approximately 
10% in the considered region of the susceptibility x- The values of the coupling constants are not very sensitive to a 
change in the continuum threshold, which gives an uncertainty of approximately 7% to the final values. 

Taking into account the uncertainties in \, Sq , and a q , the predicted values for NNa, AAcr, EScr and EEcr coupling 
constants in terms of quark-cr coupling constant read : 

9NNa/9g = 3-9 ± l.0,9AAa/9q = 1.9 ± 0.5 , <?EE<r/< = 0.4 ± 0.1 , 5sS(T / < g" q = 3.8 ± 1.0 . (38) 

To determine the coupling constants, one next has to assume some value for the quark-cr coupling constant g°. 
Adopting the value g a q — 3.7 as estimated from the sigma model we obtain 

9nN(t = 14.4 ± 3.7 , g\A<7 = 7.0 ± 1.9 , g S s<, = 1.5 ± 0.4, .g SSff = 14.1 ± 3.7 . (39) 



1 We refer the reader to Ref. l.'Ul for a detailed analysis of NNa coupling constant in QCDSR. 
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TABLE I: The scalar meson-baryon coupling constants in the SU(3) limit where the qq picture for the scalar mesons with the 
ideal mixing is assumed. 



M 


NNM 


AAM 


HEM 


EEM 


ASM 


SiVM 


ANM 


a 


14.6 


6.2 


1.3 


13.3 








fo 





-12.0 


-18.9 


-1.8 








do 


12.0 




1.3 


13.3 


6.2 






K 












-1.3 


-14.7 



Note that, the coupling constants in Eq. (fTTT^f) are denned at t = 0, i.e. gsBa = 9BBa(t = 0). As stressed above, 
the value of the susceptibility xg is less certain than the value of x- If we let XG change in a wider range, say 
6 GeV -1 < — xg < 14 GeV -1 , this brings an additional 15% uncertainty to the values of AAct and ££ct coupling 
constants but the SSct coupling constant remains intact because the sum rule is independent of the susceptibilities as 
stressed above. In order to keep consistent with the analysis in Ref. |2^, nere we also have taken Aqcd = 0.1 GeV. 
A change in the value of this parameter, say an increase to Aqcd — 0.2 GeV, does not have any considerable effect 
on SSct coupling constant, but the NNa and ££ct coupling constants are increased by approximately 8%, while the 
increase in the value of A Act coupling constant is by 5%. 

Our next concern is to investigate the SU (3) relations for the scalar meson-baryon interactions and see if the 
coupling constants above as obtained from QCDSR are consistent with these relations. The values of three coupling 
constants as determined from QCDSR together with the first equation in Eq. JSJ are sufficient to determine the three 
parameters of flavor SU(3) structure of scalar meson-baryon couplings; namely g%, g, and a s . For this purpose, 
we calculate the coupling constants in Eq. (|39[1 with the average values of the parameters; x = Xg — — 10 GeV -1 , 
a q = 0.51 GeV 2 and Sq — 2.3 GeV 2 where we obtain: 

9NNa/9g = 4.0, 5aa<t / 9q = 1.7 , g~~*/g° = 0.3 , 5sEff /g° = 3.6 . (40) 

We first assume qq structure with the ideal mixing angle 6 S ~ 35.3°, and use gNNai g~z<j an< 4 5sso- in Eq- I40[) . The 
F/(F + D) ratio, a s , can directly be calculated via the relation, 

—2a 

(g~~.v - 9NNa)/{g^a ~ gNNa) = ~ 7T~ • (41) 

1 - 2a s 

With straightforward algebra, the values of the F/ (F + D) ratio, and the octet and the singlet couplings for the qq 
picture are determined as, 

a s = F/(F + D) = 0.55 , g/fi = g N Na /g a q = 3.3 , gx/g a q = 3.2 . (42) 

Inserting a s and g into the SU (3) relations in Eq. Q and using the mixing scheme for the singlet and the octet 
couplings as in Eq. © with the value of <?i in Eq. (|42|l . we observe that the coupling constants as determined from 
QCDSR in Eq. (|40() are consistent with the SU(3) relations. This also gives gNNf — with the second equation in 
Eq. JSJl, which is justified by the non-strange content of the nucleon and by the ideal mixing scheme. In Table [I] we 
give all the scalar meson-baryon coupling constants, obtained from these relations, assuming g^ = 3.7. 

In the case of the q 2 q 2 picture, we have a quite distinct ideal mixing scheme for the scalar mesons from that for the 
qq picture. The ideal mixing angle corresponds to 9 S ~ —54.8° in this case. In this picture, we assume that the u- or 
the d-quark in the baryon couples to the qq component in the scalar meson but not to the ss component and the (qq) 
condensates are modified in the same way. Accordingly, the s-quark only couples to the ss component. Applying the 
same procedure as in the qq case with S = —54.8°, we obtain the values of F/(F + D) ratio, and the octet and the 
singlet couplings for the q 2 q 2 picture as: 

a s = F/(F + D) = 0.55 , g/g° = g N Na /g q = -2.3 , 5l /< = 4.6 . (43) 

Inserting these values into the SU (3) relations in Eq. 0) , we observe that the scalar meson-baryon coupling constants 
as found from QCDSR are consistent with the SU (3) relations, as in the qq picture. 

In Table ITT1 we present the scalar meson-baryon coupling constants g±, gs and a s in the q 2 q 2 picture, assuming 
g q = 3.7. Comparing what we have found for the scalar meson-baryon couplings in the two pictures, the value of the 
F/(F + D) ratio remains intact, as apparent from Eq. I|41|l . however, the values of g and g\ in the two pictures are 
quite different from each other. On the other hand, the BBfo couplings differ very much with regard to the structure 
of the scalar mesons. In the q 2 q 2 picture for the scalar mesons, in contrary to the qq picture, /o strongly couples to 
the nucleon due to uu and dd components it has. The strengths of the I — I couplings in the two pictures differ by a 
factor of ^/2. 
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TABLE II: Same as Table[I]but for the q 2 q 2 picture for the scalar mesons. 



M 


NNM 


AAM 


EEM 


EEA'/ 


AEM 


EAM 


ANM 


a 


14.6 


6.2 


1.3 


13.3 








fo 


10.3 


16.2 


19.7 


11.3 








ao 


-8.5 




-0.9 


-9.4 


-4.4 






K 












0.9 


10.3 



B. Sigma Meson-Baryon Coupling Constants with SU(3) Breaking Effects 

Now we turn to the effect of SU(3) breaking, where we allow m s = 0.15 GeV and / = —0.2, keeping m q = 0. We 
also restore the physical values for the parameters of baryons [32, |47fl : 

Ma = 1.1 GeV , Ms = 1.3 GeV, M E = 1.2GeV, (44) 
A A = 3.3 GeV 6 , A| = 4.6 GeV 6 , A| = 3.3 GeV 6 , 
4 = 3.1 ± 0.3 GeV 2 , Sq = 3.6 ± 0.4 GeV 2 , Sq = 3.2 ± 0.3 GeV 2 . 

The corresponding Borel windows are chosen as: 

for A , 1.0 GeV 2 < M 2 < 1.4 GeV 2 , (45) 
for 5 , 1.5 GeV 2 < M 2 < 1.9 GeV 2 , 
for E , 1.2 GeV 2 < M 2 < 1.6 GeV 2 , 

In Figs. 16181 we present the Borel mass dependence of the LHS and the RHS of the sum rules for AAct, EEct and 
SSct, respectively, with the SU(3) breaking effects. We plot each term on the RHS individually as we did in the 
SU(3) limit. We observe that the responses of the continuum thresholds for all the couplings are quite sizable. The 
contributions of the single pole terms to the SSct and EEct sum rules are large and opposite in sign to the third terms 
on the RHS. Therefore the contributions of these two terms tend to cancel each other which leads to a very stable 
sum rule for these couplings. The contribution of the single pole term for the A Act coupling is very small. 

Taking into account the uncertainties in x, Sq , and a q , the predicted values for AAct, SSct and EEct coupling 
constants in terms of quark-cr coupling constant with the SU (3) breaking effects read: 

5AAa/< = 2.0±0.5, g~~<,/g a q =0.5 ±0.1, <?WSg = 5.7 ± 1.4. (46) 

Adopting again the value g a q = 3.7 as estimated from the sigma model |46j| . we obtain 

g AAa = 7.4 ±1.9, <7HEa = 1.9±0.4, 5sSct = 21.1 ± 5.2 . (47) 

A few remarks are in order now. Comparing the values obtained from the sum rules in the SU (3) symmetric 
limit and the ones beyond the S'J7(3)-limit, we observe that the introduction of the SU(3) breaking effects does not 
change the AAct coupling constant, while the SSct and EEct couplings are modified by approximately %30 and %50, 
respectively. We also note that, the obtained value of AAct coupling constant is small as compared to the NNtt and 
NNa coupling constants. Since the a exchange gives the dominant contribution in the AA system, this suggests that 
the AA interaction is weak, in accordance with the recent experimental result. 



IV. DISCUSSION AND CONCLUSIONS 



We have calculated the AAct, SSct and EEct coupling constants which play significant roles in OBE models of YN 
and YY interactions, employing the external field QCDSR method. The coupling constants can be determined in 
terms of quark-CT coupling constant in this method. In order to compare our results with the others in the literature 
and keep as model-independent as possible, we find it useful to give the ratios of the coupling constants in the SU (3) 
limit for the average values of the vacuum parameters, 

^^ = o.43, ^1 = 0.08, ^5^=0.91. (48) 

gNNcr gNNcr gNNa 

We observe that the SSct coupling constant is more than one order small as compared to JVTVct coupling constant. 
The EEct coupling constant is at the same order with NNa coupling constant and twice as large as the AAct coupling 
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constant. We have shown that these coupling constants as determined from QCDSR satisfy the SU(3) relations which 
lead to a determination of the F/(F + D) ratio for the scalar octet. Although the scalar meson-baryon coupling 
constants depend on the picture assumed for the structure of the scalar mesons (qq or q 2 q 2 ), the F/(F + D) ratio 
remains intact in the two pictures. We would like to also note that the third terms on the RHS's of the sum rules 
which represent the responses of the continuum thresholds, affect only the values of the individual coupling constants, 
which receive contribution by the same factor. Therefore the ratios of the coupling constants and the value of a s 
remain unchanged if this term is omitted. 

All the Nijmegen soft-core OBE potential models have S > 30°, which points to almost ideal mixing angles for 
the scalar qq states. Since the ideal mixing angle for the scalar mesons has been assumed in our QCDSR calculations 
as well, it is convenient to compare our results with the ones from NSC potential models. Our result for the ratio 
9aa<t/9nn<t is half of the value found in Ref. however, it qualitatively agrees with the one from NSC89 3, 
which is QAAa / 9 n n a = 0.58. The value we have obtained for the AAct coupling constant is small as compared to 
NNa coupling constant and this implies that AA interaction is weak, since the sigma exchange gives the dominant 
contribution to this interaction in terms of OBE models. The value of the F/(F + D) ratio, which is 0.55 as obtained 
from QCDSR is about half of the values given in NSCa-f 0, which is F/(F + D) ~ 1.1. 

In order to estimate the ST/ (3) breaking in the couplings, we have restored the physical values of the parameters like 
the strange quark mass and the physical baryon masses. We observe that the SU (3) breaking effects do not change 
the AAcr coupling, while the SSc and ££<t couplings are modified largely. It is also desirable to derive the sum rules 
for the BBf and BBa couplings in order to estimate the SU(3) breaking in these coupling constants. 
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FIG. 1: The Borel mass dependence of LHS and the fitted RHS of the sum rule for AAa coupling in Eq. 1251 for Sq = 2.3 GeV 2 
and xg = X — — 10 GeV" 1 . We also present the terms on the RHS individually. Note that the LHS curve (solid) overlies the 
RHS curve (dashed). 
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FIG. 2: Same as Fig. but for the sum rule for z,Ea coupling in Eq. 13011 . 
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FIG. 3: Same as Fig.^but for the sum rule for EEcr coupling in Eq. 1311 . 
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FIG. 4: The dependence of gAAa/gq on the susceptibility \ f° r three different values of Sq = 2.0, 2.3, and 2.5 GeV 2 ; here we 
take x = XG- 
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FIG. 5: Same as Fig. 2]but for the sum rule for the gssv/gq coupling constant. 
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FIG. 6: The Borel mass dependence of LHS and the fitted RHS of the sum rule for AA<r coupling in Eq. 125H with the SU(3) 
breaking effects and for xg = X = — 10 GeV -1 . We also present the terms on the RHS individually. Note that the LHS curve 
(solid) overlies the RHS curve (dashed). 
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FIG. 7: Same as Fig. |5]but for the sum rule for HHct coupling in Eq. I3LW . 
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FIG. 8: Same as Fig.|S|but for the sum rule for SEa coupling in Eq. 13H . 



